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SUMMARY 


Investigation  of  the  turbulent -boundary- layer  flow  over  a  flat 
plate  in  compressible  flow  is  carried  out  on  the  basis  of  the  scheme 
established  in  NACA  TR  25*4-2.  By  averaging  the  Navier-Stokes  equations, 
differential  equations  for  the  mean  flow  are  obtained.  A  temperature- 
velocity  relation  follows  without  a  specified  form  of  the  length  scale. 
To  derive  the  velocity  distribution  in  the  boundary  layer,  a  choice  of 
the  length  scale  has  to  be  made.  The  temperature- velocity  relation 
reduces  to  Reynolds1  analogy  and  the  velocity  distribution  goes  back  to 
Von  Karman’s  logarithmic  law  for  the  special  case  of  incompressible 
flow. 


There  are  essentially  three  universal  constants,  arising  out  of  the 
correlations  in  the  energy  equation,  to  be  determined  by  comparing  with 
suitable  experiments  of  the  temperature-velocity  relation  at  any  known 
Mach  number  and  heat  transfer  at  wall.  The  behavior  at  other  Mach 
numbers  and  heat-transfer  conditions  may  then  be  readily  predicted. 
Because  of  the  lack  of  accurate  experimental  data,  attempts  to  carry 
out  such  determinations  are  not  included  in  the  present  report. 


INTRODUCTION 


As  part  of  an  investigation  of  Von  Karmen1 s  similarity  theory  and 
its  extension  to  compressible  flows,  the  theory  for  incompressible  flows 
was  examined  critically  in  reference  1  by  using  modern  concepts.  It  was 
found  that  the  original  form  of  the  theory  is  supported  by  modern  con¬ 
cepts.  In  reference  2,  the  theory  was  extended  to  the  case  of  compres¬ 
sible  flows  and  it  was  found  that  the  analysis  could  be  carried  through 
with  additional  approximations  but  without  any  modification  of  the  basic 
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concepts.  In  the  present  paper  the  problem  of  compressible  flow  in  a 
turbulent  boundary  layer  is  considered  on  the  basis  of  the  scheme 
established  in  reference  2.  It  is  shown  how  the  influence  of  the  Mach 
number  (e.g. ,  on  the  velocity  and  temperature  distribution)  can  be  pre¬ 
dicted  from  the  theory  after  the  constant  coefficients  in  the  theory 
are  determined  by  one  set  of  experimental  measurements. 

This  investigation  was  conducted  at  the  Massachusetts  Institute  of 
Technology  under  the  sponsorship  and  with  the  financial  assistance  of 
the  National  Advisory  Committee  for  Aeronautics. 


VELOCITY-TEMPERATURE  RELATION 


Attempts  will  now  be  made  to  deduce  relationships  among  the  mean 
quantities  and  their  distributions  within  the  boundary  layer.  To  do  so 
one  returns  back  to  the  complete  equations  of  motion  (equations  (26) 
to  (29)  of  reference  2)  instead  of  working  with  the  "localized"  equa¬ 
tions  where  the  observer  rides  with  the  mean  velocity.  These  equations 
are  averaged  so  that  only  mean  quantities  and  the  correlations  of  the 
fluctuating  quantities  appear. 

A  strictly  parallel  flow  will  be  considered.  The  justification  of 
its  application  to  the  case  of  a  boundary  layer  will  be  discussed  later. 
By  using  the  continuity  equation,  it  is  easy  to  show  that 


(1) 


when  the  mean  flow  depends  on  y  only.  Here  the  bar  denotes  the  average 
of  the  quantity  under  it  and  F  stands  for  any  function.  (See  appendix  A 
for  definitions  of  symbols.)  Thus,  equation  (26)  of  reference  2  leads  to 


dy 

_6_ 

Sy 
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2 
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There  follows  immediately,  by  integration, 

puv  =  Constant  =  tq 

PV2  =  p  -  pj  V  (3) 

pvw  =  Constant  =  T-_ 

U  /j 

where  tq  is  the  shearing  stress  at  the  plate  in  the  x-direction, 

Toz,  the  shearing  stress  at  the  plate  in  the  z-direction,  and  pj,  the 

pressure  outside  the  boundary  layer.  One  expands  the  averaged  quan¬ 
tities  in  equations  (3)  and  finds  that  triple  correlations  are  involved 
as  well  as  the  mean-flow  variables  and  double  correlations.  In  usual 
cases  the  triple  correlations  are  taken  to  be  of  higher  order  than  the 
double  correlations.  For  example. 


A  word  of  caution,  however,  may  be  in  order  here  to  warn  against  the 
overconfidence  in  such  estimations.  It  is  nothing  more  than  a  plausible 
guess  of  the  individual  terms,  and  sometimes  a  combination  of  the  indi¬ 
vidual  terms  may  very  well  invalidate  the  result.  The  discussion  in 
reference  2  of  the  dilatation  e  serves  as  a  good  example. 

Nevertheless,  it  is  perhaps  admissible  in  the  present  case  to  drop 
triple  correlations  involving  p'  on  the  basis  of  the  estimations  like 
equation  ( 4) .  If  one  concedes  that  the  effect  of  the  triple  correla¬ 
tions  is  not  negligible  but  relatively  small,  then  the  validity  of  the 
theory,  being  an  approximate  one,  might  not  be  seriously  impaired 
because  of  the  omission.  By  leaving  out  the  triple  correlations,  equa¬ 
tions  (3)  are  expanded  into  the  following: 


(p*  )20(u’v’ ),  «  p  u'v' 
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p  u'v'  +  u  p'v'  =  T Q  (5) 

P  =  Pi  +  P  (v’ )2 

+  rMl2  L  (6) 

pi  u2  y 

P  V’w'  =  Toz  (7) 

where  subscript  1  denotes  the  free-stream  quantity.  From  the  equation 
of  continuity,  one  has 


hence 


p'v'  =  Constant  =  m  (8) 

# 

where  m  is  evidently  the  mass  transfer  through  the  wall.  Likewise  from 
the  equation  of  energy  (equation  (28)  of  reference  2) 


d  — -s;  8p  — 


which  yields,  by  omitting  the  triple  correlation  p'v'T*  and  inte¬ 
grating  between  0  and  y, 


(p  vr  +  t  p'v)  -  r(p  -  Po)  = 


r>y 


Os* 


pe  dy  +  q  +  q^ 


(9) 


where  pQ  is  the  pressure  at  the  wall  (cf.  equation  (6)),  6*,  the 

effective  thickness  of  the  laminar  sublayer,  q,  the  heat  transfer  at 
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8*  _ 

the  wall,  and  q*  =1  ne  dy,  the  heat  generated  within  the  laminar 

Jo 

sublayer  (cf.  appendix  B  for  discussion  of  q*). 

One  may  now  make  the  following  observations:  In  equation  (8),  the 
turbulent  mass  transfer  m  vanishes  if  there  is  no  addition  or  sub¬ 
traction  of  mass  through  the  plate.  For  a  nonporous  flat  plate,  then, 
equation  (5)  leads  to 


p  u’v'  =  tq  (10) 

One  may  note  that  equation  (10)  is  also  the  basic  equation  in  refer¬ 
ences  3  and  4,  where  the  role  of  p'v'  was  not  mentioned.  In  refer¬ 
ence  5,  the  term  p'v'  is  kept  and  given  no  physical  interpretation. 

With  the  previous  assumption  of  small  turbulence  level 

(v')  «  1,  when  the  free- stream  Mach  number  is  of  order  unity, 

equation  (6)  reduces  to 


P-P!  (li) 

as  usually  accepted  for  boundary- layer  flows,  and  the  term  p  -  drops 

out  in  equation  (9).  It  is  seen,  however,  that  equation  (ll)  ceases  to 
be  true  if 


M,2  xtf-o(l) 

u2 


Equation  (7)  is  merely  the  statement  of  constancy  of  the  transverse 
shear  in  the  z-direction.  It  does  not  influence  the  mean  motion  in  the 
xy-plane. 

The  next  step  is  to  introduce  the  similarity  theory  and  represent 
the  fluctuations  in  terms  of  the  scales.  One  obtains  from  equation  (10) 


(12) 
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by  absorbing  the  correlation  constant  u !v 1  into  vQ,  the  similarity 

scale  for  velocity  fluctuations.  Similarly  from  equation  (9)>  with  T1 
broken  into  Tq1  and  T2*  according  to  reference  2, 


cpP 


0 


■’T1* 


aiPv0 


2  £y 

to 


+ 


q  + 


q* 


(13) 


where  aq  is  a  correlation  constant  for  the  dissipation  and  use  has 
been  made  of  the  relation 


P 


~  iW p  «  v02/t0 

Introducing  scales  0q  and  @2  for  Tq1  and  T2'  >  respectively,  and 
absorbing  the  correlations  into  proportional  constants,  09  and  ccq 

CpP^’V  +  V*T2’)  =  PVq^Zo  0  +  a3voJ 

Using  equation  (12)  and  the  definition  of  tQ  in  equations  (32)  of 
reference  2, 


—  2  cLy 

aipv°  tj  = 


v 

alTo 

us* 


du 


=  a1T0(u  -  u*) 


with  u*  representing  the  velocity  at  the  edge  of  the  laminar  sublayer. 
Hence  equation  (13)  becomes 


pv( 


{^2lo^+  a3vo2)  =  alTo(u 


u*)  +  q  +  q* 
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or 

“2  H  +  oc^Vq  =  ax(u  -  u*)  +  (l4) 

where  again  the  relations  (12)  and  equations  (32)  of  reference  2  have 
been  used. 

The  quantity  v0  may  be  put  in '  terms  of  T: 


with  p  =  ~zr.  Also,  since  q*  =  tqu*  (see  appendix  B),  one  may 
introduce  an  "effective"  velocity 


(16) 


Then  the  final  form  of  the  differential  equation  between  T  and  u  is 
reduced  from  equation  (l4)  into 


dT  pr  _ 

a2  7j=r-  +  a3P)|T  =  axue  +  q  (17) 

Equation  (17)  is  derived  without  a  specific  assumption  on  the  length 
scale  Z0.  The  constant  represents  the  heating  due  to  dissipation, 

as  a  part  of  the  external  work  tqu;  a2  represents  both  the  mixing 
phenomenon  and  the  compressibility;  and  reflects  the  combined  effect 

of  the  compressibility  and  the  dissipation.  On  a  closer  look  at  expres¬ 
sion  (l6),  it  is  revealed  that  the  factor  (l  -  j  arising  from  the 

heating  effect  is  originated  from  the  idealized  situation  of  an  entirely 
turbulent  outer  boundary  layer  with  similarity  and  an  inner  laminar 
viscous  sublayer.  In  the  sublayer  all  the  energy  supplied  by  the 
external  work  must  go  into  heat  through  dissipation,  hence  the  constant 
unity  appears.  In  the  turbulent  layer,  part  of  the  energy  is  trans¬ 
ferred  or  turned  into  turbulent  energy,  therefore  only  a  portion  a^_ 

is  dissipated.  If  a  transition  region  between  the  two  idealized  regions 


8 


NACA  TN  25^3 


had  been  assumed,  there  would  be  no  discontinuity  in  the  picture  but 
a  term  of  the  order  ^1  -  a^Ju*  would  still  be  present  in  the  heating 
because  of  dissipation. 1 

Equation  (17)  niay  be  nondimensionalized  by  writing 


6  =  t/t-l 

a  =  *el* 1 

cf  =  To/|  Pi 

cq  =  q  /  P1  ul(To  "  Tl)cv 


(18) 


where  subscript  1  again  refers  to  the  free- stream  value  and  TQ  is 
the  temperature  at  the  wall.  Then  there  follows 


^■3  =  M-^A-jU  +  A2^)+  Q  (19) 

du  '  ' 


■4?o  see  this,  one  could  assume  a  linear  transition  of  the  dissi¬ 
pation  parameter  such  that 


lie  =  td  0  <  y  ^  cx6*,  cx  <  1 

'^*1/ 


1  + 


(al  “  X) 


u  -  u. 


uCo  -  UC 


du 

dy 


,  c-,5*  ^  y  ^  Cq  >  1 


du  •  „  ^ 

=  a*|Tn  — ,  CoS*  ^  y 

±  u  dy 

where  c^b*  <  y  ^  CpS*  is  the  transition  region  and  uc^  and  uc  are 
the  respective  velocities  at  the  ends  of  the  transition.  Then  the  term 


corresponding  to  ( 1  -  — \u*  is  (l  -  — 

v  ai/  v  ab 


u_  +  u~„ 
C1  c2 
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*  where 
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Ai  =  ?R 


^1 

a2 


Q  = 


(20) 


Equation  ( 19)  is  the  fundamental  equation  of  the  present  theory  for  the 
variation  of  mean  temperature  with  mean  velocity.  In  the  spirit  of  the 
local  similarity,  since  the  turbulence  pattern  is  obtained  by  letting 
the  observer  ride  with  the  mean  local  velocity,  it  seems  quite  clear 
that  the  Mach  number  effect  at  most  would  be  an  indirect  one.  The 
constants  ap,  a^,  and  a 3,  being  essentially  correlation  constants, 

ought  to  be  nearly  independent  of  the  Mach  number.  The  effect  of  Mach 
number  on  the  mean  temperature -velocity  relationship  is  largely  in  the 
first  term  on  the  right-hand  side  of  equation  (19),  where  Mp2  stands 

before  the  parenthesis.  One  may  also  expect  the  constants  A2  and  Q 
to  vary  in  some  way  with  the  Mach  number  through  the  frictional  coef¬ 
ficient  cf  and  the  heat-transfer  coefficient  cq.  However,  the  varia¬ 
tions  of  cp  and  Cq  with  Mach  number  are  known  to  be  rather  slow. 

For  an  approximate  theory,  it  might  be  sufficient  to  regard  Ap,  Ap, 

and  Q  also  as  universal  constants  independent  of  the  free-stream  Mach 
number.  Likewise,  the  constants  Ap  and  A2  might  further  be  taken 

to  be  independent  of  the  heat-transfer  situation  at  the  wall. 

There  are  consequently  three  parameters  ap,  02,  and  a3  to  be 

introduced  in  the  problem  besides  the  boundary  conditions  of  the  mean 
flow  at  the  wall  and  in  the  free  stream.  These  parameters  are  related 
to  the  turbulence  mechanism  and  must  be  empirically  determined  by  com¬ 
paring  with  experimental  data.  Once  determined  from  suitable  experi¬ 
ments  at  low  speeds  with  known  heat  transfer  at  the  wall,  they  enable 
one  to  predict  the  behavior  at  any  other  Mach  number  or  heat-transfer 
condition,  provided,  of  course,  that  the  turbulence  level  is  still 
small  and  that  the  Mach  number  is  not  excessive. 

In  the  following,  'equation  (19)  will  be  integrated  first  for  the 
case  of  subsonic  flow  and  then  for  the  general  case.  For  the  subsonic 
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case  the  solution  can  he  expressed  in  an  ascending  series  of  M^2,  and 
a  few  terms  would  he  sufficient  to  compare  with  low-speed  tests  to 
deduce  the  values  of  the  universal  constants.  The  integral  for  the 
general  case  is  useful  mainly  for  the  prediction  of  supersonic  boundary 
layer  after  the  universal  constants  are  determined. 

Case  of  subsonic  flow.-  Assume  the  solution  to  be  expanded  in 
ascending  powers  of  M^2: 

6  =  9(0)  =  M129(1)  +  M^9(2)  +  .  .  .  +  M2n9(n)  +  .  .  .  (21) 


hence, 


W=  e(0)1^2  +  Mlf!  0(O)"1/20(lJ]  + 

M14U  0(O)'l/^][®(2)  -  i  e(O)'l0(l)^|  + 

Mi6ie(or1/2]f(3)  -  ie(orle(i)e(2) 


|0(of20(i)3]  +  • 


(22) 


Substituting  into  equation  (19)  and  equating  powers  of 
a  set  of  differential  equations  for  the  functions  ^(o) 
Thus,  for 

M-l0: 

de<°)  ,  „ 
dft 

M!2: 

^v~-2W/£ 

M^: 

i2)  '  l^e(o)-l/2e(i) 

M-l6: 

i3’  “  1  A2e(0)-l/2[f(2) 

1  n  - 1, 

'4e(o) 

2n 


(23) 
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and  so  forth.  For  boundary  conditions  one  has  for  the  incompressible 
case  vith  no  heat  transfer  0  =  0(q)  =1  when  ft  =  u^,  where 


This  result  holds  for  any  Mach  number;  therefore, 


0(l)  =  0(2)  =  •  •  •  =  0  when  u  =  u-^ 

With  these  conditions  the  integrals  of  equation  (23)  are  found  to  be 


and  so  forth.  When  the  wall  is  insulated,  Q  =  0,  and  the  functions 
0(1)  y  0(2)  >  •  •  •  are  simply  polynomials  Of  u.  For  convenient 
reference,  a  few  of  them  are  listed  as  follows: 


and  so  forth. 
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Aq  interesting  result  to  be  noted  is  that  the  solution  9(o)  in  * 

equation  (24)  evidently  coincides  with  what  is  usually  known  as  Reynolds' 
analogy,  where  a  similar  transfer  mechanism  is  assumed  for  both  the  » 

momentum  and  the  heat  transfers.  Here  the  same  results  are  obtained 
because  of  the  fact  that.  the  turbulent  exchange,  proportional  to  the 
temperature  gradient  dT/dy,  dominates  the  situation  at  very  small  Mach 
numbers.  A  similar  transfer  mechanism  follows  the  assumed  similar  tur¬ 
bulence  pattern.  The  solutions  and  so  forth  now  give  the  cor¬ 

rection  to  Reynolds'  analogy  at  higher  Mach  numbers. 

General  case.-  A  brief  outline  will  now  be  given  for  the  general 
case.  Rewrite  equation  (19): 


Jits  _  w: 

\du  2 


i J  and  u.  Using  the  standard  method,  put 
:  F(u)^A-1^1  U  +  | j  (27)  * 


which  is  now  homogeneous  in 


■ 


AiMi2 


^  Q 
u  +  — 

2 


(26) 


Then  F  is  to  he  solved  from 


F  dF 


du 


A1M12  V2  _  W2 
2  2 


1 

2  2 


(28) 


Integrating  equation  (28),  one  gets 


Al^1  |  =  C^F  -  k1)Sl(F  -  k 2)S2 


(29) 
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where 


Together*  equations  (27)  and  (29)  form  a  parametric  representation  of 
the  function  9  =  0(11).  Applying  the  condition  at  the  outer  edge* 
namely*  0  =  1  when  u  =  uj_*  one  has 


and 


The  complicated  way  in  which  the  constants  are  entangled  renders  their 
evaluation  from  empirical  data  highly  tedious.  Besides*  the  proper  form 
of  equation  (29)*  where  every  term  is  real*  depends  on  the  sign  and 
magnitude  of  the  constants  and  is  best  to  be  directly  integrated  from 
equation  (28)  once  the  constants  are  known. 


VELOCITY  AND  TEMPERATURE  DISTRIBUTIONS  IN 
BOUNDARY  LAYER 


To  determine  the  distribution  of  the  mean  velocity  and/or  the  mean 
temperature  across  the  boundary- layer  thickness*  it  is  necessary  to  bring 
in  the  definition  of  the  length  scale  l0  as  a  function  of  y.  In 

reference  2*  two  alternative  expressions  (equations  (33)  and  (39))  have 
been  given  for  the  length  scale,  laying  emphasis  on  the  momentum  and  the 
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energy  equations,  respectively.  Integrating  either  will  lead  to  a 
velocity  distribution,  from  which  a  temperature  distribution  may  be 
calculated  by  means  of  the  results  in  the  section  "Velocity-Temperature 
Relation."  The  distributions  evaluated  from  the  two  expressions  will, 
in  general,  be  different.  It  seems  that  only  by  comparison  with 
experiments  can  one  say  whether  both  lead  to  essentially  the  same 
result  or  one  form  i’s  preferable  to  the  other  in  certain  particular 
cases. 

The  two  definitions  are,  nevertheless,  of  the  same  general  form: 


I  cc 

0  <w 


(32) 


with  X  standing  for  the  variable  u  in  equation  (33)  of  reference  2 
and  for-  T  in  equation  (39)  of  reference  2.  Without  specifying  X, 
it  follows  from  equation  (32)  of  reference  2  that 


,,  _  i  du 

Vo  a  ‘'o 

dy 


=  a 


dy 


where  is  a  proportional  constant.  Rewriting, 


Hence,  by  integration  and  substitution  of 


equation  (15), 


,  dX 
10g<=  to  ■  ‘1 


t 


(33) 


fl 

0 


(3*0 
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The  temperature- velocity  relation  may  then  he  used,  yielding 


dX 

dy 


=  G(u) 


where 


G(u)  =  exp 


'fi  [V 

SJff. 


(35) 


Finally,  one  may  obtain  the  distribution  by  integrating  equation  (35) 


dX 

G(u) 


dX  d£  du 
d@  du  G(u) 

For  definition  (33)  of  reference  2,  X  =  u, 


7 

PdX  d£ 
~J  d0  du 


(36) 


r-  f-S 

J  0(5 


(37) 


for  definition  (39)  of  reference  2,  X  =  9, 


y  = 


d£  du 
du  G(u) 


(38) 


Again  one  may  first  solve  the  subsonic  case  for  a  more  explicit  expres¬ 
sion,  which  must  be  reducible  to  the  usual  logarithmic  distribution  for 
incompressible  flow.  Then  a  brief  discussion  of  the  general  case  for  anv 
Mach  number  will  be  taken  up.  < 


Case  of  subsonic  flow.-  The  function  G(u)  must  first  be  evaluated 
before  integrating  equation  (36).  By  definition  (35), 
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after  nondimensionalizing  with  equation  (8).  Using  equation  (22), 
one  has 


where 


u 


U1 


=  90(u)  +  M-^cp-^u)  +  M^cp2(u)  + 


<P0(S)  - 


°u 


e(o)'l/2  a“ 


Tl(K)  - 


cp2(^)  = 


'U 


'u 


U1 


e(0)‘3/2  0(1) 


du 


u- 


(l-  6/-,  re, 


2fl  -5/2 

Q  °(l)  0(O) 


|  e(2)e(0)" 3/2) 


cp3(u)  = 


Pu  , 

J 


e(l)e(2)S(0) 


-5/2 


I  e(3)e(o) 


-3/2 


1?  e(D39(0)'7/2 1  « 


!) 


and  so  forth,  with  the  0^Ts  defined  by  equation  (24)  with  heat 

transfer  at  wall  and  by  equation  (25)  with  an  insulated  wall.  Sub¬ 
stituting  into  equation  (39);  one  may  therefore  put 


(4o) 
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G(u)  =  g0(u)jl  +  Mx2 g1(u)  +  M^ggCu)  +  .  .  .J 


where 


g0(u)  =  exp  (ax 


tp0(u)  +  a2j 


«1(“>  * 


'P12(S) 


e3<a)  ‘  al  \f^ 


<P3(“)  *  al\J^T  %reK'2!“)  *  5^  ’’I3'21 


V  (31) 


J 


Inverting, 


1 

G(u) 


+  M12g_1(u)  +  M1^g_2(u)  + 


/ 


where 
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and  so  forth.  Equation  (42)  may  now  be  used  to  evaluate  equations  (37) 
and  ( 38 ) : 

(l)  Case  (a)  X  =  u: 


du 

G(S) 


=  Yu^0(u)  +  M12Yu^1(u)  +  M14Yu^2(u)  +  .  .  . 

where  yr  is  a  reference  station,  u^.  is  the  value  of  u 
at  yr,  and 


'u 


du 


ur  Bb(S) 


Y  m 


■u,l  - 


a  Si(u) 


g0(u) 


du 


gi2(S)  -  gg(g) 
e0(u 


It  will  be  of  interest  to  apply  equation  (43)  to  the  incompressible 
case,  namely  M-^  — >  0.  The  distribution  is  then 
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For  simplicity,  one  may  take  yr  =  8,  8  being  the  thickness  of  the 

? 

boundary  layer,  so  that 

"r  “  U1 

u  -  U!  =  jL  -  1 
U1 

.  With  given  by  equation  (24), 


^Actually  this  procedure  does  not  give  the  best  fit  to  experimental 


velocity  distributions,  since  the  similarity  concept  is  not  valid  near 
the  outer  edge  of  the  boundary  layer.  Instead  of  5,  the  length  scale 
should  be  taken  as  proportional  to  v/uT,  where  uT  is  the  frictional 
velocity.  For  general  discussions,  however,  there  ought  to  be  little 
dif f erenc  e . 
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where 


The  constant  a£  may  be  determined  from  the  usual  condition  at  the 
wall:^ 


as 


The  other  constant  a^  has  to  be  evaluated  by  matching  with  experi¬ 
mental  distribution.  From  the  definition  (33),  it  corresponds  to 
Von  Karman's  universal  constant. 


One  may  further  develop  the  solution  (44)  into  a  power  series  in 
Q  for  cases  where  the  heat  transfer  is  relatively  small.  After  some 
manipulation,  the  final  result  is 


where 


^As  a  possible  refinement,  the  laminar  sublayer  may  be  assumed  to 
have  a  linear  velocity  profile,  and,  instead  of  infinity,  the  velocity 
gradient  at  the  wall  may  be  prescribed  as 


y 


0 


However,  this  step  complicates  the  practical  calculation  and  must  be 
justified  by  experience. 
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Applying  the  condition  of  infinite  velocity  gradient  at  the  wall/  one 
finds  easily 

b2  =  -6  (46) 

For  the  case  of  an  insulated  wall,  there  follows 


u 


(47) 


which  is  precisely  a  form  of  the  logarithmic  law,  matching  with  the 
experimental  distribution  at  the  assumed  edge  of  the  boundary  layer. 
Such  a  form  has  been  suggested  by  Dryden  (reference  6)  in  1935,  where 
the  frictional  velocity  Uy  and  Von  Harman' s  constant  K  were  used 

instead  of  a-]_  and  cf  in  equation  (47).  By  recalling  the  definition 
one  verifies  readily  that  ' 


Equation  (47)  then  reduces  exactly  to  Dryden’ s  form 

u  -  uj  =  -  ^  loge  |  (49) 

One  may  also  conclude  that  the  logarithmic  law,  suitably  modified, 
indeed  can  be  approximately  true  even  with  the  presence  of  a  very  small 
amount  of  heat  transfer,  since  equation  (45)  shows  that  the  first-order 
correction  for  Q  =  0  merely  amounts  to  replacing  the  factor  in  equa¬ 
tion  (49)  by  a  slightly  different -one.  If  the  heat  transfer  is  appreci¬ 
able,  equation  (45)  predicts  that  the  logarithmic  law  in  the  form  of 
equation  (49)  would  break  down.  There  would  be  considerable  interest 
if  equation  (45)  could  be  tested  for  different  values  of  Q. 

It  is  apparent  that  each  of  the  fractions  Yu  in  equation  (43) 

may  be  developed  into  a  power  series  in  Q  as  above.  In  general, 
therefore,  the  right-hand  side  of  equation  (46),  which  gives  the  velocity 
distribution,  is  in  the  form  of  a  double  series  in  the  parameters  Q 

and  ,  involving  constants  a-^  and  a2  which  are  to  be  determined 
once  and  for  all,  in  fact,  by  comparing  with  the  incompressible 
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insulated  case  as  above,  both  a^_  and  ag  are  given  in  terms  of 
B-|_  and  B2« 

I 

For  its  practical  importance,  the  results  for  the  insulated  case 
at  subsonic  Mach  numbers  are  explicitly  given  as  follows:  Equations  (4o) 
are  integrated  by  using  equations  (25), 


and 

bn  = 


cp0(u)  =  u  -  ux 

cp^(u)  =  bQ  +  b-^u  +  b2U^  +  b^u^ 

/  ~  \  iv  /v>Q  "5  <viii  /v>\) 

q>2(u)  =  Cq  +  c^u  +  CgU^  +  c^u-J  +  c^u  +  c^u 

so  forth,  where 
At  ^  o  Ao  ^  p 

!  ■  T  U1  '  T  U1 


(50) 


,  -^1  ~  2  -^2  ~ 
^1  =  X  Ul  +  “  Ul 


bo  = 


.tSL 


"bn  —  - 


12 


c0  = 


Cn  — 


"  ££  Al2^l5  -  ^  AlA2Sl4  -  in  -  (jsA^  +  I  A22)^ 

^  Al2^!4  +  h  A1A2S13  +  J  A22S12  +  (J  A1A2  +  |  A22)*l  ^  ^ 


c2 

c3 


3  .  ,  ~  2  3,2~  1.,  1  . 

-  -  £  ^i-^2ul  ~  a  a2  U1  ”  ~7  ^1^2  ~  £  Ap 


I  a  2 


8  *  A  16 


"  TZ  A12u12  "  |  A1A2u1  +  ^  A22 


16 

1_ 

16 


8 


C1*  -  S  A1A2 

c5  -  -i-  Ai2 
5  160  1 
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and  so  forth.  The  functions  g  of  equations  ( 4l)  become 


6q(u)  _  exp  Jb-^u  -  u^J  +  a^j 


g^S)  =  B1cp1(u) 


g2(u)  =  B1j^»2(a)  +  |  B1cp12(u)~| 


and  so  forth.  Hence , 


Yu,0  =  B2 


Yu,l  "  B2 


Yu,2  =  B2 


B2  -  exp  [-Bl(u  -  uj]  j 

B2  jjBl(u  "  ^10  +  dUu  +  dl2u2  +  d13^)- 

(d10  +  dliei  +  d12^12  +  d13^13)j- 

b2  |jexp  jjBl(S  "  al)  "  a2j(d20  +  d21^  +  d22’^2  + 

^23U^  d24u  ^  +  {^20  +  ^21ul  + 

d22%2  +  d23%3  +  deV3flJ*'  +  d25^15  +  d26^16j| 


and  so  forth,  where 


d  =b  f  h  ■  21°2 ,  6b3 

10  0  B1  Bl2  Bl3 

du  =  bx  +  p-  +  % 

11  1  Bx  Bl2 

d12  =  b2  +  57- 


d13  =  B3 
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and 


&2,  n  = 


6-n  .  .  , 

=  V—  (m  +  n)l  hm 

*  n!  -n  m* 

m=u  Bj 


The  constants  hm  are  defined  as  follows: 


Bi  o 

ho  =  -c0  +  —  bQ 


hl  =  "cl  +  BlVl 


h2  =  "Cg  + 


^(b-,2  +  2h0b2) 

ll(blb2  +  b0u3 > 


b3  c3  +  Bllblb2  +  ^0^3) 


h4  =  -ci,.  +  ^(b22  +  2bib3) 

=  -c^  +  ^1^2^ 3 

h  _  B1  -u  2 

h6  '  ~  b3 

With  these  evaluations,  the  distribution  for  the  insulated  subsonic  case 
becomes  finally, 


(55) 


(56) 


where 


loge  |  =  -Bl^=  -  *  M12f1(u)  +  U1kf2Cn) 

fi(s)  -  -  Y2  -  v) 


+  #  *  • 


f. 


(S)  -  E  (42n  -  I  T“  -  V) 

£0  \  l+m^n  /'  7 


(57) 
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so  that  a  modification  of  the  definition  of  B2  leads  to  the  same 
distribution  (45).  Thus  the  two  choices  of  X  are  equivalent  to  each 
other  even  for  the  distribution  in  incompressible  flow. 

For  the  compressible  case  with  an  insulated  wall,  the  right-hand 
side  of  equation  (60)  starts  with  Y Now 


* 


by  taking  yr  =  5  and  using  equations  (23)  and  (25).  Integrating, 


Ye,i  '  "b2  (a2  +  by  +  Al^Je 


-Bl(u-ui) 


~  \2  +  %  +  AlUl 


Substituting  into  expression  (59) >  one  has 


y  -  5  =  +  o(m/) 

The  condition  of  infinite  velocity  gradient  at  y  =  0  again 
determines  Bg, 

S  =  -B2U±2(a2  +  ^  4.  Alux 
Hence  the  distribution  may  be  written  as 
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Equation  (62)  indicates  that  in  the  absence  of  heat  transfer  at  the 
wall,  the  choice  of  X  =  6  leads  to  an  additional  term  to  the  usual 

logarithmic  law  even  when  «  1  (hut  not  zero).  The  additional 

term  is  of  more  weight  nearer  the  wall  where  u^uj  differs  appreciably 
from  unity.  Noticing  further  that  A3JA2  is  a  controlling  parameter, 

fl  =  al/a3 

a2  ’  FT2 

I 

by  equation  (20),  one  may  say  that  the  predicted  deviation  from,  the 
usual  logarithmic  law  depends  on  the  relation  between  the  heat  gen¬ 
erated  by  dissipation  and  the  turbulent  transfer  of  such  heat  (and 
due  to  compressibility).  Ordinary  theories  assume  the  dissipation  to 
be  negligible,  equivalent  to  putting  A-j_  =  0.  Then  the  additional 

term  vanishes. 

The  apparent  contradiction  of  equations  (45)  and  (62)  when  only 
the  lowest -order  terms  in  both  are  retained  is  resolved  by  the  obser¬ 
vation  that,  in  complete  form,  one  should  have 


y  -  8  =  QYU,0  +  M1%,1  +  •  *  •  (^3) 


so  that  the  lowest-order  terms  should  be  taken  according  to  whether 
Q  «  M-j_2  or  the  converse  is  true.  If  both  are  of  the  same  order, 

neither  should  be  left  out  and  the  final  expression  certainly  will 
contain  an  additional  term  analogous  to  that  in  equation  (62). 
Physically,  it  is  obvious  that  if  a  great  amount  of  heat  conduction 
is  present  at  the  wall,  it  certainly  would  overshadow  the  distributed 
heating  due  to  dissipation  (and  compressibility)  at  lower  Mach  numbers. 
On  the  other  hand,  when  there  is  no  external  heating,  the  dissipation 
becomes  the  only  predominant  factor. 

The  Mach  number  effect  in  equation  (62)  can  be  worked  out  without 
difficulty,  but,  involving  no  new  features,  the  details  are  omitted  in 
this  report. 

General  case.-  The  general  case  will  be  only  briefly  treated,  as 
the  integration  depends  on  the  values  of  the  constants  A^,  A2,  and  Q. 
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(Cf.  equations  (28)  to  (30).)  Still  starting  with  equations  (36) 
and  (39)>  one  substitutes  ^  by  equation  (27) 


where 


G(u)  =  exp 


nu 


B- 


du 


+  a- 


F(fc) 


rAnM 


rVii 


tt  + 


=  exp  -B 


d F 


+  a. 


F  A1M1  F2  .  A^M1  F 


-  1 


=  B 


F  -  ki^l 

3\f  -  k^  ) 


Pi  = 


A1M12  k2  "  kl 


B3  =  e 


a2 


F(ui)  -  k2 


F(ui)  -  k]_ 


Pi 


(64) 


Hence ,  for  X  =  u, 


y  - 


yr  = 


du 


«r 


G(S) 


Pf 


B 


du 


'f  -  kh'Pl 


Fr  ^ 


-  k- 


dF 


(65) 
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with  du/dF  given  by  equations  (28)  and  (29).  Similarly,  for  X  «  6, 


PF 

1_  d©  du 

Bo  dil  dF 

3  jFr 


/F  -  kiVPl 

\F  -  k 2 ) 


dF 


(66) 


with  d0/du  to  he  evaluated  from  equations  (27)  and  (29).  The  con¬ 
stant  of  integration  B3  is  determined  from  the  condition: 


as 


The  constant  has  previously  been  identified  with  u^K. 


CORRELATION  OF  THEORY  WITH  EXPERIMENTS  AND  DISCUSSIONS 
Determination  of  Arbitrary  Constants 


There  are  two  sets  of  arbitrary  constants.  In  the  temperature- 
velocity  relation,  one  has  three  constants  aq_,  and  representing 

essentially  combinations  of  the  correlations,  which  are  assumed  uni¬ 
versal  with  regard  to  the  free- stream  Mach  number  and  the  heat-conduction 
conditions  at  wall.  In  addition  there  is  the  mean  velocity  u*  at  the 
edge  of  the  laminar  sublayer  separating  the  turbulent  boundary  layer 
from  the  wall.  The  second  set  arises  out  of  the  integration  for  the 
mean  velocity  distribution  within  the  boundary  layer,  consisting  of 
B1  and  B2  in  the  section  "Velocity  and  Temperature  Distributions  in 

Boundary  Layer."  The  external  conditions  defining  the  mean  flow. .-are 
the  following:  The  free-stream  velocity,  pressure,  temperature,  the 
friction  at  wall,  and  the  temperature  at  wall  (or  the  amount  of  heat 
conduction  from  the  wall). 

^-th  given  external  conditions,  the  constants  can  be  determined 
by  comparing  with  accurate  experimental  results.  The  significant  point 
is  that,  once  determined,  a  universal  law  for  all  Mach  numbers  and  heat 
conduction  is  established.  No  other  ad  hoc  assumptions  need  to  be 
introduced  for  special  cases. 
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It  is  then  only  necessary  to  conduct  a  few  low-speed  experiments 
with  prescribed  heat  conduction  at  the  wall  and  measure  ^the  mean  tem¬ 
perature  and  velocities  at  points  across  the  section.  Elias1  data 
(reference  7)  were  in  this  category  but  unfortunately  were  not  pre¬ 
sented  in  enough  detail  to  suit  the  present  purpose.  One  may  naturally 
also  deduce  from  high-speed  measurements  with  somewhat  more  work. 

Wilson  (reference  8)  and  Ladenburg  and  Bershader  (reference  9)  published 
the  results  for  a  supersonic  free  stream.  Unfortunately,  again,  their 
data  were  not  sufficient,  because  only  one  of  the  two  variables  u 
and  T  was  measured  in  either  case.  Both  assumed  the  isoenergetic 
relation  to  hold  for  evaluation  of  the  other  variable. 

There  are  arguments  for  the  validity  of  the  isoenergetic  relation, 
either  from  the  practical  reason  that  the  effective  Prandtl  number  is 
nearly  unity  or  from  some  assumption  on  the  mixing  lengths  for  momentum 
and  energy,  such  as  Ferrari1 s  (reference  5).  Even  admitting  the  argu¬ 
ments,  no  isoenergetic  law  would  follow  when  the  dissipation  term  is 
kept  in  the  energy  equation.  In  reference  2,  it  has  been  shown  that 
the  dissipation  terms  should  be  retained  for  a  consistent  theory,  and 
one  could  only  regard  the  success,  if  any,  to  be  essentially  an 
empirical  one  for  the  Mach  number  or  heat  conduction  involved.  Such 
empirical  results  might  be  useful,  however,  for  the  approximate  evalua¬ 
tion  of  the  constants  a^,  and  0,3,  at  least  as  a  guidance  to  the 

orders  of  magnitude. 

The  constant  u*  is  a  more  troublesome  one.  In  appendix  B,  it 
is  shown  to  be  dependent  on  the  Mach  number  as  well  as  another  empirical 
constant  defining  the  extent  of  the  turbulent  layer.  The  correct  determi 
nation,  therefore,  is  very  tedious  and  perhaps  not  warranted  because 
of  the  approximate  nature  of  the  theory.  A  possible  way  is  to  deter¬ 
mine  its  value  from  incompressible  flow  and  neglect  the  dependence  on 
Mach  number.  The  error  introduced  can  be  assessed  only  by  comparing 
with  experimental  results . 

The  determination  of  B^_  and  B2  is  done  by  matching  with  the 

empirical  distribution  of  the  mean  velocity  or  temperature.  In  the 
section  "Velocity  and  Temperature  Distributions  in  Boundary  Layer"  it 
is  found  that  B2  is  essentially  determined  by  the  condition  that 

du/dy — >  00  as  y  — ►  0,  and  B^  corresponds  to  u r|K  in  the  usual 

incompressible  case  (cf.  equation  (47)  and  (48)),  K  being  Von  Karman1  s 
universal  constant.  Approximately,  within  the  range  of  moderate  Mach 
numbers,  B^  might  be  taken  to  be  independent  of  Mach  number  without 

serious  error. 


NACA  TN  2543 


31 


Some  Experimental  Results  of  Turbulent 
Boundary  Layer 

A  Brief  examination  of  the  velocity  distribution  in  the  incompres¬ 
sible  case  may  be  made  to  check  with  the  logarithmic  law  predicted  by 
similarity  theory.  A  quite  common  form  for  turbulent -boundary-layer  flow 
has  been  the  l/7-power  law.  Recently,  there  are  the  careful  measurements 
in  the  National  Bureau  of  Standards^  and  some  others  such  as  those  by 
Hama  (reference  10 ).  Typical  examples  from  these  works  at  various 
Reynolds  numbers  have  been  plotted  on  semilogarithmic  paper  as  fig¬ 
ure  1.  It  is  found  that  outside  the  "viscous"  layer  near  the  wall, 
all  the  curves  tend  to  show  a  more  or  less  straight  portion,  of  approxi¬ 
mately  the  same  slope,  the  extent  of  the  straight  portion  decreasing 
as  the  Reynolds  number  decreases. Such  a  phenomenon  is  in  accordance 
with  the  expectation  on  the  basis  of  the  similarity  theory,  because  a 
smaller  Reynolds  number  brings  in  the  viscous  effects  to  a  greater  degree, 
causing  a  merging  of  the  higher-  and  lower-frequency  parts  of  the  energy 
spectrum.  The  slope  was  predicted  to  be  proportional  to  uT  and  hence 
should  change  only  slowly  with  Reynolds  number,  again  verified  by  the 
experiments.  For  larger  values  of  y/S,  the  velocity  curve  deviates  from 
the  straight  portion  on  the  semilogarithmic  plot  and  may  be  approximated 
by  a  power  law.  In  this  part  of  the  boundary  layer,  the  turbulence  is 
nearly  isotropic  and  therefore  the  mean  flow  is  not  expected  to  follow 
closely  the  logarithmic  relation.  The  empirical  l/7-power  law  is  also 
included  and  shows  good  agreement  on  the  whole. 

For  the  compressible  case,  the  works  of  Wilson  (reference  8)  and 
Ladenburg  and  Bershader  (reference  9)  have  been  quoted  before.  The 
velocity  distributions  given  are  also  plotted  in  figure  1  for  comparison. 
It  should  be  noted  that  in  deducing  the  velocity  distribution  from  their 
measurements,  both  made  use  of  the  isoenergetic  relation.  Wilson’s 
measurements  were  done  with  a  pitot  tube  and  registered  actually  the 
Mach  number  variation  within  the  boundary  layer,  computed  from  the 
stagnation  pressures.  A  discussion  was  made  on  the  error  introduced 
by  the  assumed  isoenergetic  relation,  with  the  conclusion  that  the 
maximum  error  amounts  to  about  3  percent  for  Mj  =2.0.  Now  Ladenburg 
and  Bershader  measured  by  means  of  an  interferometer  and,  therefore. 


^As  yet  unpublished.  The  authors  are  obliged  to  Dr.  Schubauer, 
chief  of  the  Aerodynamics  Section  in  NBS,  for  furnishing  the  data. 

^Hama's  results  show  more  waviness  than  those  of  NBS,  partially 
because  his  Reynolds  number  is  quite  low,  and  are  omitted  in  figure  1. 
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have  recorded  actually  the  temperature  variation.  By  coincidence  the 
two  sets  of  experiments  have  some  of  the  tests  made  under  almost 
identical  conditions,  like  the  following: 

Reference  9'-  =  143°  C,  M^  =  2.3>  Rg  ~  10^ 

Reference  8:^  T-^  =  150°  C,  M-^  =  2.03,  Rg  ~  leP 

After  correcting  for  the  Mach  number  difference,  one  may  estimate  the 
temperature  distribution  in  Wilson's  case  from  Ladenburg  and  Bershader's 
data  and  evaluate  the  velocity  from  the  measured  local  Mach  number.  In 
correcting  for  the  Mach  number  difference,  assume  a  modified  isoenergetic 
law  to  hold, 


~  »  1  +  aMx  1 1 
T,  V 


a  being  a  constant,  so  that 


Then,  letting  subscripts  refer  to  the  respective  experiments, 


(T\  ~ 

/T\ 

( M-i  -  f  M-i 

\  1  /  W  \  ^  /  L  and  B 

tai  ~ 

\Ti/l  and  B 

\Ti/l  and  B 

(V) 

'  1  'L  and  B 

and,  finally, 


^Table  I,  station  6,  of  reference  8. 
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The  following  results  are  then  obtained: 


y/8 

o.io4 

0.217 

0.321 

0.477 

0.686 

("K)w 

.665 

.751 

.820 

•  903 

• 

VO 

03 

PO 

1.31 

1.24 

1.20 

1.12 

1.05 

(u|%)w 

.762 

.836 

.899 

.958 

1.003 

wNiL 

.769 

.838 

.888 

.943 

.990 

where  the  subscripts  WW  refer  to  the  value  given  by  Wilson  using  the 
isoenergetic  law.  It  is  interesting  to  see  that,  at  least  for  this 
particular  case,  the  isoenergetic  law  does  give  a  very  close  approxi¬ 
mation  as  checked  by  the  independent  temperature  measurement.  However, 
whether  the  same  agreement  would  result  in  other  cases  cannot  be  readily 
ascertained. 


Variation  of  Skin  Friction  with  Reynolds 
and  Mach  Numbers 


In  the  present  theory  the  skin  friction  was  introduced  as  an 
external  condition  such  as  would  influence  the  turbulence  pattern 
through  its  effects  on  the  mean  distributions.  One  may  stretch  the 
theory  a  little  to  yield  an  equation  of  the  skin  friction  involving  the 
chordwise  Reynolds  number  and  the  Mach  number,  as  Prandtl  and  Von  Kdrm&n 
had  done  for  the  incompressible  case  (references  11  and  12).  A  momentum 
thickness  for  the  entire  boundary  layer  may  be  derived,  which  is 
related  to  the  skin  friction  by 


T 


—  —2  dfl 

P1U1  s 


The  viscosity  is  introduced  to  form  a  length  scale  with  the  frictional 
velocity  uT  and  a  chordwise  Reynolds  number  appears.  When  the  univer¬ 
sal  constants  in  the  present  theory  are  known,  the  skin-friction  relation 
thus  derived  will  also  give  explicit  Mach  number  effects.  A  point  to  be 
noted  is  that  the  velocity  distribution  should  refer  to  the  frictional 
velocity  uT,  instead  of  u£  at  the  outer  edge  y  =  5.  For  simplicity 

it  might  be  necessary  to  make  approximations  but  there  seems  to  be  no 
great  difficulty. 


Massachusetts  Institute  of  Technology 

Cambridge,  Mass.,  December  27,  1950* 
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APPENDIX  A 

SYMBOIS 


1’  2 


constants,  defined  by  equations  (33)  and  (39); 
respectively 


^0)loV  *  *  * 
c0>cl*  *  *  * 


coefficients  in  expansions  for  cp^  and  cp2 
defined  in  equations  (51) 


C-p,C 


f  >  q 


'P'  v 


coefficients  of  skin  friction  and  heat  transfer; 
respectively;  defined  in  equations  (l8) 

specific  heats  at  constant  pressure  and  constant 
volume ;  respect ively 


d10>dll>  * 
d20>d21>  * 


fl,f2 


g0'gl> 


E>q  ,  ^  y  • 


^0* "  * 
k 

kx,k2 

l0 


coefficients  in  expansions  for  Yu^p  and  Yu?2j 
defined  hy  equations  (54)  and  (55),  respectively 


dilatation 


(du  bv_  8w\ 

8x  8y  8z  / 


correction  functions  for  Mach  number  effect  on 
logarithmic  lav;  defined  hy  equations  (57) 

functions  in  expansion  for  G;  defined  hy 
equations  (4l) 

functions  in  expansion  for  l/G;  defined  hy 
equations  (42) 

constants  defined  hy  equations  (56) 

coefficient  of  heat  conductivity  of  fluid 

constants  defined  hy  equations  (30) 

similarity  scale  of  length 


m 


rate  of  mass  transfer  at  vail 


0 


0 


r 


m 
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PjPl 

q 

* 

q 

s1js2 

t 

to 

U,Y,V 


u* 


u 


mean  pressure  and  mean  pressure  in  free  stream, 
respectively 

rate  of  heat  transfer  at  wall 

rate  of  amount  of  heat  generated  in  laminar  sublayer 

constants  defined  by  equations  ( 30) 

time 

similarity  scale  of  time  for  fluctuations 

velocity  components  in  x-,  y-,  and  z-directions, 
respectively 

similarity  scale  of  velocity  fluctuations 
modified  mean  velocity,  defined  by  equation  ( l6) 

magnitude  of  mean  velocity  at  edge  of  laminar 
sublayer 

nondimensionalized  modified  mean  velocity  ^u  =  u~|u^) 


U1 


value  of  ft  evaluated  at  outer,  edge  of  turbulent 
boundary  layer 


yr 


value  of  u  evaluated  at  a  reference  station  within 
boundary  layer 

frictional  velocity 


Cartesian  coordinates;  x,  axis  in  direction  of  plate 
and  free  stream;  y,  axis  normal  to  plate;  and  z, 
axis  parallel  to  leading  edge  of  plate 

reference  station  in  boundary  layer 

constants  defined  by  equations  (20) 


Bi,B2 


constants  defined  by  equations  (45) 


constant  defined  by  equations  (3l) 
function  defined  by  equation  (27) 

function  defined  by  equation  (35) 
Mach  number  of  free  stream 
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Q 

R 

*5 

T 

% 

U 

X 


Yu,0'Yu,l'  *  • 

Ye,o’Ye,i’  •  * 
al,Ct2,a3 

3 

Pi 

7 

5 

5* 

€ 

0 

0l'e2 

%reur 
d(2)>  *  '  * 


constant  defined  by  equations  (20) 

gas  constant  in  equation  of  state 

Reynolds  number  based  on  boundary- layer  thickness 

temperature 

mean  temperature  at  wall 

mean  velocity  component  in  x-direction 

function  defining  similarity  scale  of  length,  as 
in  equation  (32) 

distribution  functions  for  mean  velocity,  defined 
by  equations  (43) 

distribution  functions  for  mean  velocity,  defined 
by  equations  (60) 

constants  representing  correlations  in  averaged 
equation  of  energy 

constant  ^R  ^2.  ^ 

constant  defined  by  equations  (64) 

ratio  of  specific  heats 
thickness  of  boundary  layer 
thickness  of  laminar  sublayer 
rate  of  dissipation 

nondimensionalized  mean  temperature,-  defined  by 
equations  (l8) 

similarity  scales  of  temperature  fluctuation 

functions  in  expansion  for  0,  defined  by 
equations  (24)  or  (25) 


t3 

v 

P 

T 


T0>  roz 


VV  •  •  • 


momentum  thickness  of  boundary  layer 
coefficient  of  viscosity 
coefficient  of  kinematic  viscosity 
density  of  fluid 
shearing  stress 

shearing  stress  at  the  vail  in  direction  of  main 
stream  and  parallel  to  leading  edge,  respectively 

functions  in  expansion  for  F,  defined  by 
equations  (4o) 


The  subscript  1  denotes  quantities  in  the  free  stream.  Barred 
quantities  always  represent  mean  values;  primed  quantities  represent 
fluctuations. 
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APPENDIX  B 

DISSIPATION  IN  THE  VISCOUS  SUBLAYER 


Imagine  the  Boundary  layer  to  be  idealized  to  consist  of  two  parts: 
The  major  part  is  entirely  turbulent,  extending  from  the  outer  edge  to 
a  certain  distance  close  to  the  wall.  From  there  on  the  viscous  effects 
take  over  completely,  and  a  viscous  sublayer  is  formed.  In  the  viscous 
sublayer,  the  flow  is  laminar  and  free  of  turbulence. 

Let  the  thickness  of  the  viscous  sublayer  be  8*.  The  energy 
equation  within  the  sublayer  then  degenerates  into 


(Bl) 


since  all  the  correction  terms  vanish  for  the  laminar  flow  with 
u  =  u(y),  v  =  0.  The  first  term  is  the  viscous  dissipation;  the 

second  is  the  heat  conduction.  Noting  that  h  =  To  =  Constant,  one 

dy 

gets  by  integration  between  0  and  8*, 


—  ,  dT 

Tnu*  -  k  — 
0  dy 


y=0 


*1 


(B2) 


assuming  conduction  to  be  comparatively  negligible  at  y  =  8*.  Obvi¬ 
ously  t0u*  is  the  heat  generated  in  the  viscous  sublayer  by  dissipa¬ 
tion,  and  -k  —  is  the  heat  conducted  into  the  fluid  from  the 

dy|y=o 

wall.  In  the  notation  of  equation  (9) ,  they  are  the  quantities  q* 
and  q,  respectively.  Therefore,  q^  is  the  resultant  heat  exchange 

between  the  viscous  and  the  turbulent  layers. 

It  is  important  to  estimate,  at  least  approximately,  the  magni¬ 
tude  u*  for  the  present  theory  (cf.  equation  (l6)).  To  do  so  the 
first  step  is  to  estimate  the  thickness  8*.  Usually  the  thickness 
of  the  sublayer  is  expressed  by 
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6*uT 

-  =  h,  a  constant  (B3) 


In  the  actual  case  there  is  a  transition  layer  where  the  viscous  and 
turbulent  shears  are  of  the  same  order  of  magnitude.  The  value  of  h 
in  the  incompressible  case  varies  from  approximately  10  to  30, 
depending  on  whether  the  sublayer  is  taken  to  be  strictly  laminar  and 
has,  consequently,  a  linear  velocity  profile,  or  the  sublayer  is  to 
consist  of  the  entire  thickness  where  noticeable  departure  from  the 
turbulent  velocity  profile  occurs.  In  the  idealization  adopted  in  this 
section,  the  transition  layer  is  omitted  and  the  boundary  of  the 
"viscous  layer"  should  therefore  lie  somewhere  within  the  actual  tran¬ 
sition  layer.  By  such  reasoning,  it  seems  that  an  average  value  of  h 
of  the  order  of  20  is  probably  adequate  for  the  purpose. 

It  may  be  seen  that,  with  such  a  concept,  the  value  of  h  ought 
to  vary  with  the  Mach  number.  For,  one  may  regard  the  sublayer  as  owing 
its  existence  to  the  fact  that  the  viscous  shear  becomes  an  appreciable 
part  of  the  total  shear.  Suppose,  then,  that  h  is  defined  to  satisfy 


du 


dy 


y=8* 


Toa 


(b4) 


a  being  a  constant  less  than  unity.  The  quantity  du/dy  may  be 
estimated  by  using  the  turbulent  velocity  distribution.  For  instance, 
in  the  subsonic  case  if  one  uses  equation  (57), 


u  =  ^  loge  y  +  B  +  M-^Cu)  +  M-^f^u)  + 
K 


dy  KyL  11  J 


Hence,  by  substituting  into  equation  (b4). 


h  -  —(l  +  +  .  .  .) 

ocK  \  / 


4o 
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where  U*  =  —  u*  (cf.  equation  (l6)).  If  h^  is  the  value  for  the 
al 

incompressible  case,  there  follows 


—  =  1  +  M-^f/tu*)  +  .  .  .  (B5) 

hi 


Since  the  right-hand  side  involves  u*,  successive  approximations  might 
he  necessary  for  obtaining  the  solution. 

Having  determined  5*,  one  may  find  u*  by  again  applying  the 
turbulent  distribution. 
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